Fl proton structure function from the unified DGLAP/BFKL approach 



<3\ 

o 
o 



oo 



43 

i ' 

<d : 

43 . 



(N 

m 
o 

On 
O 



X 



K. Golec-Biernat 1 ' 2 -* and A. M. Stasto 3 ' 4 ' 1 ^ 

'institute of Nuclear Physics PAN, Krakow, Poland 
s Institute of Physics, Rzeszow University, Rzeszow, Poland 
3 Penn State University, Physics Department, University Park, PA 16802, USA 
4 RIKEN Center, Brookhaven National Laboratory, Upton, NY 11913, USA 

(Dated: May 8, 2009) 

We compute the longitudinal proton structure function Fl from the kr factorization scheme, 
using the unified DGLAP/BFKL resummation approach at small x for the unintegrated gluon 
density. The differences between the kr factorization, collinear factorization and dipole approaches 
are analyzed and discussed. The comparisons with the HERA data are made and predictions for 
the proposed Large Hadron-Electron Collider are also provided. 
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I. INTRODUCTION 

The longitudinal nucleon structure function Fl, mea- 
sured in the deep inelastic lepton-nucleon scattering, is 
proportional to the cross section for the interaction of 
the longitudinally polarized virtual photon with a nu- 
cleon. This observable is of particular interest since it is 
directly sensitive to the nucleon gluon distribution. In 
the naive quark-parton model Fl vanishes (the Callan- 
Gross relation). This is due to the quark spin 1/2 and the 
fact that the struck quark has limited transverse momen- 
tum in the naive parton model. In the QCD improved 
parton model, however, the gluon interactions cause the 
average quark transverse momentum to grow with 
increasing value of the (minus) photon virtuality Q 2 . As 
a result, Fl acquires a nonzero leading twist contribution 
proportional to a s (Q 2 ). At small values of the Bjorken 
variable x, Fl is driven mainly by gluons through the 
transition g — > qq. Therefore, it can be used for the ex- 
traction of the gluon distribution in a nucleon providing 
a crucial test of the validity of perturbative QCD in this 
kinematical range. 

The experimental determination of Fl is in general 
difficult and requires a measurement of the inelastic cross 
section at the same values of x and Q 2 but for different 
center-of-mass energy of the incoming beams. This was 
achieved at the DESY electron-proton collider HERA by 
changing the proton beam energy with the lepton beam 
energy fixed. The structure function Fl was measured 
both by the HI [1] and ZEUS [2] collaborations in the Q 2 
range of 12 — 90 and 24 — 110 GeV 2 , respectively. 

At small x, the nucleon structure functions receive 
large logarithmic corrections coming from resummation 
of large powers of a s lnl/x. This procedure goes be- 
yond the standard collinear factorization and is achieved 
by the use of the kr factorization formalism [3, 4] 
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with the unintegrated gluon density found as a solution 
to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [5-7] or 
Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution 
equations [8-11]. Since the small x expansion receives 
large corrections at higher orders, resummation at small 
x is in general necessary in order to obtain predictions 
which are in agreement with data. 

The objective of this paper is the calculation of Fl 
within the fcr factorization formalism using the unin- 
tegrated gluon density obtained from the Kwiccihski- 
Martin-Stasto (KMS) approach [12], which provides 
a convenient framework for the unification of the 
conventional Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) [13-15] and small x BFKL evolution equa- 
tions. From the point of view of the small x hierarchy, 
the KMS approach includes important effects of higher 
order resummation. In addition, we systematically ana- 
lyze the relation between this approach and the collinear 
and dipole approaches, investigating the role of differ- 
ent contributions to Fl in various kinematical regions. 
We especially emphasize the role of the exact gluon kine- 
matics in the kx factorization formulae and demonstrate 
numerically that this kinematics have a sizable effects on 
the predictions for Fl, and thus, on the extracted gluon 
density. We compare our computations with the exper- 
imental data at small x from the HI [1] and ZEUS [2] 
collaborations. 

The paper is organized as follows. In Sec. II we re- 
call the fcr factorization formalism for the longitudi- 
nal structure function. In Sec. Ill we review the uni- 
fied BFKL/DGLAP approach for the unintegrated gluon 
density which includes important resummation effects at 
small x. In Sec. IV we discuss the relation of the kr fac- 
torization with the collinear and dipole approaches. In 
Sec. V we present a systematic numerical analysis of the 
various approaches and compare them with the HERA 
data, as well as provide the extrapolations to the LHcC 
(Large Hadron-Electron Collider [16]) energies. In the 
last section we summarize our conclusions. 
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FIG. 1: Quark box diagrams in the photon-gluon fusion process used in the fey-factorization formula. 



II. F L FROM THE k T FACTORIZATION 
APPROACH. 

In the limit of the high center-of-mass energy, or equiv- 
alently at small values of Bjorken x, the nucleon structure 
functions can be computed from the kr factorization ap- 
proach [3, 4]. Basic diagrams for the boson-gluon fusion 



which are taken into account in the high energy limit are 
depicted in Fig. 1. The gluon is off-shell with its virtu- 
ality dominated by the transverse momentum kx = k, 
see Fig 1. Since we are interested in Fl, the photon 
is longitudinally polarized. Thus, in the kr factoriza- 
tion approach the longitudinal structure function is then 
given by 



9 



dV a s (/j, 2 )p 2 (1 



^ 2 \D lq D 2q 



f(-,k' 

\ z 



(1) 



where the denominators D\ q and D 2q read 



D lq = n 2 + f3(l~f3)Q 2 + 
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the quark transverse momentum kt = k, the shifted 
transverse momentum is given by tt' = n — (1 — 0)k, 
and the argument of the unintegrated gluon density 
f(x/z, k 2 ) is defined to be 



x _ 
z 



X 1 



k' 2 + m 2 




The variable (3 is the corresponding Sudakov parameter 
appearing in the quark momentum decomposition, 

(5) 

K 2 \ 

(6) 




(l-P)Q 2 J 
with the following light-like base vectors 



M 2 x 



P =P 



q = q + xp 



(7) 



where M denotes the nucleon mass, m q is the quark mass 
which we keep nonzero only for charm quark, p is the 



target-proton four-momentum and q is the virtual photon 
four-momentum. The argument in the strong coupling 
constant is taken to be 



2 a 
fl = K 



(8) 



The integration over the gluon virtuality k 2 in eq. (1) 
needs special care in the low momenta region, k 2 < k 2 ~ 
1 GeV 2 . We will discuss this important element of our 
presentation in the forthcoming sections. 

The function f(y,k 2 ) is the unintegrated gluon dis- 
tribution, which in the small x limit is related to the 
conventional (integrated) gluon distribution g(y, /i 2 ) by 



dk 2 



(9) 



The integration limits in (1) are constrained by the con- 
dition x g < 1 while the condition x g > x is automat- 
ically satisfied from eq. (4). We note that in the strict 
high-energy limit the argument of the unintegrated gluon 
distribution would be set to the Bjorken x. This is also 
the usual procedure in the dipole picture approach which 
we discuss in Sec. IV. Here, we take into account the ef- 
fects of exact kinematics which results in the shift of the 
gluon x g to larger values than x. This is related to the 
fact that the energy needed to produce the qq pair is 
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non-negligible even when the total center-of-mass energy 
is very large. Although this effect is non-leading in the 
leading logarithmic small x approximation, it is never- 
theless numerically quite important, as we will illustrate 
in Sec. V. 

It also has been shown in the dipole picture that by in- 
cluding the exact kinematics in the argument of the gluon 
distribution, the transverse size of the quark-antiquark 
dipole is no longer conserved [17, 18], see Sec. IV B. 

III. UNIFIED DGLAP AND BFKL EQUATIONS 

The main input to the kr factorization formula is the 
unintegrated gluon distribution f(y, k 2 ). At small x, this 
distribution can be found from the solution to the BFKL 
[5-7] or the CCFM equations [8-11]. These equations 
give predictions for the unintegrated gluon density as a 
function of the transverse momentum squared k 2 and x 
(and also an external scale Q in the case of the CCFM 
equation) provided x <C 1. A more rigorous approach 
which includes the operator definitions of unintegrated 
gluon densities is presented in [19]. Here, we will use the 
unintegrated density obtained from the solution to the 
set of unified BFKL and DGLAP equations, which in- 
cludes the small x resummation effects. The full formal- 
ism, called the KMS approach, was constructed in [12] 
and here we only review main elements of this approach. 

A. Equation for the unintegrated gluon density 

In the KMS approach [12] one constructs the evolution 
equation for the unintegrated gluon distribution function 



which includes the leading order BFKL kernel with the 
kinematical constraint and the DGLAP part of the split- 
ting function. The solution to the leading order BFKL 
equation is well known [20], giving very fast growth of 
the gluon density with the decreasing value of x. It can 
be recast into the symbolic form: 

f(x,k 2 )~x-\ (10) 

where A ~ 0.5 is the so-called hard Pomeron intercept. 
This fast growth was shown to be incompatible with the 
experimental data which exhibit an effective intercept 
A ~ 0.3. 

The next-to-leading corrections to the BFKL equation 
turned out to be very large [21-23] and it became imme- 
diately apparent that additional resummations are nec- 
essary. In the KMS approach one uses the leading log- 
arithmic approximation for the BFKL kernel but with 
substantial modifications. One of them is the kinemati- 
cal constraint [24, 25] which accounts for a large portion 
of the next-to-leading order corrections. It is important 
to stress that this constraint goes beyond the next-to- 
lcading order in logarithms of x, and is responsible for 
partial resummation of the small x series [26] . Another 
modification in the KMS approach is a non-singular part 
in 1/z of the splitting function P gg in the BFKL kernel 
in addition to the already included singular part. In a 
series of papers [27-29] (see also [30-35]) it was shown 
in detail how these modifications generate higher order 
terms in the small x expansion. 

The final equation in the KMS approach, which takes 
into account all the modifications mentioned above, has 
the following form 




where the strong coupling constant a s = a s N c /ir. The 
first term on the r.h.s is the non-perturbative input to 
be specified below, while the second term contains the 
leading logarithmic BFKL kernel with the kinematical 
constraint given by the theta function. The third term 
contains the DGLAP splitting function P gg (z) without 
the singular term in z, and the last term is the contribu- 
tion from the quark to gluon transition with X being the 
singlet quark distribution. The input function for this 



integral equation is chosen to be 

/ (0) (*,* a ) = ^ j\zP gg (z)- z g(- z ,k 2 ). (12) 

Note that the special form of the input is dictated by the 
fact that eq. (11) only involves f(x,k 2 ) in the pertur- 
bative domain, k 2 > fcg, where Kq is a non-perturbative 
cutoff taken to be equal k 2 = 1 GeV 2 . The gluon in- 
put (12) is provided by the conventional gluon distri- 
bution xg{x, Kq). This guarantees consistency with the 
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DGLAP evolution equations since the input in both ap- 
proaches is exactly of the same. In this way, the necessity 
of parametrizing the unintegrated gluon distribution in 
the non-perturbative regime, k 2 < k 2 , is avoided. 

B. Equation for the singlet quark density 

In the KMS approach the equation for the unintegrated 
gluon density was supplemented by the second equation 
for the quark density. These two equations formed the 
coupled system of equations (similarly to the DGLAP 



where the quantities k', D\ q , £>2 9 and x g were defined in 
Sec. II. 

The singlet quark momentum distribution contains 
both the sea and valence quarks. The contribution to the 
singlet quark distribution was calculated differently de- 
pending on the region of the transverse momenta. There 
are three regions of interest for k 2 and k' 2 : 

1. non-perturbative: k 2 ,n' 2 < k$ 

2. strongly ordered with low gluon transverse mo- 
menta: k 2 < fcg < k' 2 

3. perturbative: k 2 > k 2 . 

In the non-perturbative region, the sea contribution is 
assumed to be dominated by the soft-Pomeron exchange. 
This part is parametrized phenomenologically in the fol- 
lowing form 

S^ix) = S£ + Sf + S? (15) 

with the soft pomeron contribution 

SZ = SS = 2S? = Cp*-°- 0S (1-*) 8 . (16) 

The second contribution comes from the region of small 
transverse momenta of the gluon, k 2 < k 2 < k' 2 . In this 



equations) for the functions f(x,k 2 ) and T,(x,k 2 ). The 
sea quark contribution was evaluated from the kr factor- 
ization theorem at small x, 

Sf\x,Q 2 ) = J'^J ^-S^(z,k 2 ,Q 2 )f( Xg ,k 2 ) 

(13) 

where S^ 30 ^ describes the quark box contributions, 
shown in Fig. 1, which implicitly includes the integration 
over the quark transverse momentum k. The explicit 
expression for (13) reads 



(14) 



I 

region the strongly ordered approximation for the quark- 
gluon transition is applied and the relevant contribution 
is given by the following formula 

S^\x, Q 2 ) = f *L S^\z, k 2 = 0, Q 2 ) U X -, kf) 

Jx z Z \ Z / 

(17) 

where the on-shell approximation, k 2 = 0, is applied to 
evaluate S^°*\ 

In the perturbative domain, k 2 > k 2 , the quark con- 
tribution is evaluated from the fcx factorization formula. 
The final expression for the singlet quark distribution is 
taken to be the sum of the contributions from the three 
discussed regions 

S = (stT+S^ l) + S^ + (S^ + S^) + V. (18) 

Note that for the charm evaluation we did not use the 
soft contribution since we assume that charm is generated 
dynamically from gluons and that there is no soft or non- 
perturbative charm contribution. 

Using the kr factorization and all the terms discussed 
above, one finds the the final equation for the singlet 
distribution S in the KMS approach: 



S^ k) (x,Q 2 ) =&J d -^j\ P J d 2 .'a s y 
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FIG. 2: Fl from the kr factorization approach. Dashed-dotted (black) line: contribution from the boson-gluon fusion box with 
gluon transverse momenta k > fco; dashed (black) line: non-perturbative input from the gluons via collinear formula; dotted 
(black) line: contribution from the quarks. Solid (red) line is the sum of all contributions. Gluon kinematics is exact. 



S^fc 2 ) = S^\x) + £ r^S^\zM 2 = ^r4)-g(-,ki) + V{x,e) (19) 



where 5 (soft )(o;) is given by eq. (15) and the uds subscript 
indicates that the additional S — > S term is only included 
for the light quarks. Eqs. (19) and (11) form a set of 
equations to be solved in the KMS approach. 



The final formula used for the calculation of the longi- 
tudinal structure function Fl from the kr factorization 
formalism with the KMS approach reads: 
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— 1^' 



^0(k 2 >k 2 )j\p J A'c« s «^-ffl 2 (i-i] f(x g ,k 2 



a s (Q 2 ) J 4 f 1 dy fx 
v I 3 J x y \y 



(20) 



r 



where the cutoff for the gluon momentum fc 2 = 1 GeV 2 
and the non-perturbative input was taken to be 

yg(y, k 2 ) = N(l - yf (21) 

with the parameters N = 1.57 and (3 — 2.5 which were 
found from a fit to the HERA data on the structure func- 
tion i<2, using the same approach. In the forthcoming, we 
will discuss in detail the relation of the kr factorization 
formula (20) to those from other approaches. 

It is interesting to see what is the magnitude of the sep- 
arate contributions to Fl in cq. (20). In Fig. 2 we show 
the breakdown of Fl into the contributions from gluons 
from the boson-gluon box (first term)), quarks (second 
term) and the non-perturbative gluon input (third term). 
The non-perturbative input stays nearly constant as a 
function of x and Q 2 . At low x the dominant contribution 



is from the gluon density in the kx factorization frame- 
work. However, this contribution is small at x > 0.01 and 
at small values of Q 2 . This is due to the kinematic ef- 
fects since phase space for the gluon emissions shrinks in 
this regime. On the other hand, the quark contribution 
is non-negligible in the same regime. 



IV. RELATION OF THE k T FACTORIZATION 
TO OTHER APPROACHES 



A. Relation to the collinear factorization approach 

The standard collinear factorization formula for the 
longitudinal structure function reads 



J 



Fl(x, Q 2 ) = 



a s (Q 2 ) 



dy 

y 



yg(y,Q 2 ) 



dy 

y 



F 2 (x 1 Q 2 



(22) 



r 



where quark masses in this formula are neglected. Thus, 
Fl has two contributions: originating from quarks and 
proportional to F 2 , and from gluons and proportional to 
the integrated gluon distribution g(y,Q 2 ). 

It is instructive to illustrate that the on-shell limit of 
the kr factorization formula (1) is compatible with for- 
mula (22). By the on-shell limit we mean the approxima- 
tion in which the transverse momentum of the gluon k 2 is 
much smaller than the virtuality of the photon, k 2 <§C Q 2 . 
To this aim, we start by expanding the expression under 
the integral in eq. (1), 



~ — 

2 \D lq 



1 



(23) 



in powers of k 2 /Q 2 . We retain only the leading term, 
proportional to k 2 , and drop all the higher powers of k 2 . 
Explicit expressions for the denominators D\ q and D 2q 



read 
D lq 



= k' 2 + 2(1 - 0)k' ■ k + (1 - f3) 2 k 2 
+ (3(l~f3)Q 2 + m 2 



D 2q - k' 2 -2/3k' -k + f3 2 k 2 + (3{l- 13)Q 2 + m. 
After expanding the denominators in k 2 we obtain 
2(1 -/?)«/ -ft 



1 — 1 
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D lq 


" a; 


1 — 1 
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D 2q 
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0{k 2 ) 



20k' ■ k 



0(k 2 ), 



where 



D„ 



k' 2 + 13(1 - I3)Q 2 + 



(24) 



is independent of the gluon transverse momentum k. 
Here, we need to keep only terms linear in k. Therefore, 
expression (23) to the first order in k 2 reads 



1 



1 



4cos 2 </>«/ 2 k 2 



(25) 
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FIG. 3: Comparison between the kr factorization and collinear factorization predictions for Fl as a function of Q 2 for two 
values of x: 10 -3 , (left) and 10 -4 (right). The solid lines are the kr factorization prediction with exact kinematics, eq. (20), 
while the dotted lines correspond to the dipole approximation, that is x g — > x in the gluon density. The dashed lines show the 
collinear factorization predictions, eqs. (37,38). 



The integration cPk' can be written as ^dK' 2 d<p and one 
can perform the azimuthal integration over the angle <j>. 
The dependence on k 2 is now only in the unintegrated 
gluon distribution. Since we have assumed the strong or- 
dering in the transverse momenta, we can easily perform 
this integration using the following definition 



and carefully set the integration limits. Using relation 
(29) we can write the inequality 



yg(y,^ 2 ) = / -grf(y,k 2 ) 



/3(1-/3)Q 2 (!_i)_ m 2> , 
and since 1 > > 0, we have 



(31) 



(26) 



m 2 q x 



where the scale /x 2 ~ Q 2 . Note that, formally the inte- 
gration over the k in eq. (1) is over all scales. However, 
we expand the kx factorization formula for small values 
of k 2 /Q 2 and therefore, assuming that transverse mo- 
menta are small. Using relation (26), we can rewrite the 
approximate form as 



(32) 



^From inequality (32) we obtain the lower limit for y 



y > x 1 



Am 2 



Q 2 



(33) 



F. 



(on - sholl) (x,Q 2 ) = 2^-Ve 2 



J d/3 J dn a s (Q ) \ s now convenient to make another change of variables: 



X /3 2 (l-f3) 2 —yg(y^ 2 ) (27) 
1 



where now 



^=2 +A < 
and using inequality (32), we finally obtain 



(34) 



x 1 



k' 2 + m 2 

- P)Q 2 



(28) 



since we have dropped the ratio k 2 /Q 2 in x g , see eq. (4). 
It is convenient to change the integration variables in 
eq. (27) from k' 2 to y, 



4 Q 2 (y-x) 



< A < 



4 Q 2 (y-x) 



(35) 



As a result, we obtain the following expression for the 
on-shell limit of the k factorization formula 



K ' 2 = 



[3(1 - (3)Q 2 (| - l) 



— 11— m„ 



(3(1 -(3)Q 



(29) 
(30) 



F 



(on— shell) 



(z,Q 2 )=2]T e 



/(i) _ 9 — £ /( 2 ) 



(36) 
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where 



j(d = f 1 Yx_£ 

9 ^ Jx q y \yj V 2/ 



'i- 



■yff(2/,Q 2 ) 



Q 2 (y-a:) 



(37) 



and 



^ 2) 



y \y 



x In 



1 



1 - 



Q 2 (y-x) 



A.m? n x 



yg(y,Q 2 ) (38) 



Q 2 (y-x) J 

with the lower cutoff on the integration equal to 



4m 2 

= x H 1 



2 



(39) 



Formula (36) together with eqs. (37) and (38) is the on- 
shell approximation derived from the kr factorization in 
the presence of quark masses. In this derivation we also 
assumed that the argument of the coupling constant is 
equal to the external scale /z 2 ~ Q 2 . It is straightforward 
to verify that the above expressions coincide with the 
gluonic contribution of the standard massless collinear 
formula (22) in the case when the quark masses vanish. 
Therefore, the collinear formula arises as a leading twist 
part of the kr factorization formula, and the second term 
in eq. (36) contains a part of the higher twist contribution 
proportional to 1/Q 2 . 

Let us emphasize that in order to obtain the limit of 
the kx factorization formula consistent with the collinear 
factorization, it was crucial to take the exact kinematics 
for the argument of the gluon density, x g , see formulae 
(1), (4) and (28). The numerical relevance of the above 
substitution in the kx factorization formula is shown in 
Fig. 3. In this figure we show the Q 2 dependence of Fl 
obtained within the kx factorization approach with and 
without the exact kinematics, as well as the computa- 
tion within the collinear approach. The differences be- 
tween the collinear and kx factorization approaches are 
not significant although there is distinctive difference in 
the slope of the Q 2 dependence. On the other hand, the 
differences due to the exact gluon kinematics are quite 
substantial and persist even in the small x regime. 



B. Relation to the dipole approach 

The dipole representation for the inclusive cross section 
can be computed from the fc-r factorization formula. It is 
obtained after the Fourier transformation of expression 
(1) from the space of quark transverse momenta k into 
the the space of the transverse coordinates, r. It is im- 
portant to note that one also needs to perform the small 



x approximation in the argument of the gluon density in 
formula (1), 



x . 



(40) 



This is obviously justified only in the limit of very small 
x. In this way the Fourier integrals over the k variable 
in (1) can be easily performed. 

In the case of the longitudinal structure function, 



F T , = 



Q 2 



47r 2 a er 



(41) 



where the longitudinally polarized photon-proton cross 
section ul reads 

x J^a s f(x,k 2 )(l-e-^ k )(l-e^ k ). (42) 
Here k is the gluon transverse momentum, the variable 



Q = - f3)Q 2 + m 2 



(43) 



Ko is the Bcssel function and r is the transverse size 
of the qq pair. The expression with the integral over 
k in cq. (42) is proportional to the dipole-proton cross 
section, a, which characterizes the interaction of the qq 
color dipole with the proton: 



<j(x, r) 



2tt 

Y 



^a s f(x,k 2 )(l-e-^ k )(l 



jr-k\ 



(44) 

Notice that such a Fourier transformation is only possible 
if the substitution (40) is done. Otherwise, by including 
the exact kinematics in the argument of the gluon distri- 
bution, the transverse size of the quark-antiquark dipole 
is no longer conserved [17, 18, 36]. 

In [37] a systematic analysis of the twist expansion (i.e. 
the expansion in powers of 1/Q 2 ) in the dipole model ap- 
proach was performed. Using the GBW saturation model 
[38] for the dipole-proton cross section, a complete hierar- 
chy of the twist series has been established. The analysis 
has been extended in [39] to a saturation model which in- 
cludes the DGLAP evolution [36]. The higher twist terms 
are proportional to the nonlinear terms in the gluon den- 
sity. Consequently, the leading twist part is the term 
which is linear in the gluon density. The leading twist- 2 
part in the dipole picture in the case of the cross section 
for the longitudinally polarized photon reads [37] 



(twist-2) _ 



Aita e 



a s xg(x, Q 2 ) 

Q 2 



(45) 



or for the longitudinal structure function 

F (twi s t- 2)= 1 ^eta s xg(x,Q 2 ). (46) 
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FIG. 4: Comparison between the collinear factorization formula (22) with only the gluon term and massless quarks (solid lines) 
and the dipole twist-2 predictions (dashed lines) as a function of Q 2 for two values of x: 10 -3 (left) and x = 10 -4 (right). 



As we noted above, the dipole picture can be recovered 
from the kr factorization formula (1) by taking the high 
energy limit and making the substitution (40). In order 
to recover the leading twist-2 expression in the high en- 
ergy limit, it is therefore enough to consider the gluon 
part of the collinear formula (22) and replace the inte- 
grated gluon density by 

yg(y,Q 2 )^xg(x,Q 2 ), (47) 

which allows to pull it outside the integral. The integral 
over y can be then performed exactly 

where we assume that x <C 1. Combining the above re- 
sult with eq. (22), we obtain the twist-2 contribution to 
Fl in the small x approximation given by eq. (46). There- 
fore, the leading twist term from the dipole approach is 
actually an approximation to the leading twist collinear 
factorization formula (22). In Fig. 4 we demonstrate nu- 
merically that the values of Fl obtained from eq. (46) 
are larger than those from the collinear factorization for- 
mula. 

Summarizing, we have demonstrated how the collinear 
and dipole approaches are related to the kx factorization 
formula (1). The collinear formula is obtained upon ex- 
panding the fc-r factorization formula in powers of k 2 /Q 2 
and retaining the lowest order in this expansion. The 
exact gluon kinematics has to be taken into account in 
that procedure. On the other hand, the dipole approach 
is obtained from the kr factorization expression in the 
limit when x is very small, which amounts to approxi- 
mating the gluon longitudinal momentum fraction x g by 
the Bjorken x. 



V. COMPARISON WITH THE HERA DATA 

We start our numerical analysis by comparing the cal- 
culations of the longitudinal structure function Fl per- 
formed within the collinear formalism and the kx factor- 
ization formalism with and without the exact kinematics. 
To be precise, for the collinear calculation we use the for- 
mulae (36,37,38) with the light quark masses set to zero 
and the charm quark mass equal to m c = 1.5 GeV. There 
is also quark contribution in the form of the second term 
in Eq. (22) which is proportional to F%. We emphasize 
that we have fitted only the Fi data, thus the calculations 
for Fl are then the absolute predictions. 

In Figs. 5 and 6 we show the calculations obtained us- 
ing the collinear approach and the results from the k^ 
factorization formalism with the exact kinematics. The 
Fl structure function is plotted as a function of x in 
bins of Q 2 . The data shown in Figs. 5 and Fig. 6 are 
from ZEUS [2] and HI [1] experiments respectively. The 
agreement between the experimental data and our cal- 
culations is good. In the case of the kr factorization 
with the exact kinematics the results are rather close to 
the ones obtained from the collinear approach. For the 
kinematics range of the HERA data we do not see any 
significant differences for Fl between these two differ- 
ent factorization schemes. The only regions where the 
results differ is the region of very small x and high Q 2 , 
where the kx factorization with exact kinematics tends 
to give higher values, and the region of small Q 2 (below 
10 GeV 2 ) where the factorization-based approach falls 
below collinear one, compare Fig. 3. This result is con- 
sistent with the previous observations [40] which show 
that the significant difference between the high energy 
and the collinear factorizations is more pronounced for 
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the transverse structure function [12]. 

We stress though the importance of the exact kinemat- 
ics in the evaluation of the gluon density. The collincar 
and kx factorization approaches give very similar results 
only in the case when the gluon density is evaluated at 
x g in the kx factorization formula. In Figs. 7 and 8 we 
show also the calculation where in the kx factorization 
the argument of the gluon density equals the Bjorken x. 
Clearly, the results which do not take into account the 
exact kinematics are much higher than those with the 
exact kinematics. This is understandable as we are tak- 
ing into account that finite energy has been used for the 
production of the qq pair, and as a result the argument 
of the gluon density x g > x . We see that the differences 
are quite pronounced, they are typically larger than the 
differences between the collinear and the kx factorization 
with the exact kinematics. The differences are also visi- 
ble in the plots of Fx as a function of Q 2 . It is interesting 
that the differences do not seem to vanish as a function 
of Q 2 . This difference can be of course accounted for 
by changing the gluon density. The results for the F L 
and F 2 structure functions can be made consistent within 
the two calculations (with and without exact kinematics) 
at the expense of having different normalizations for the 
gluon density. 

We have checked that the contribution from the kx 
factorization which is proportional to the gluon density 
is about 2 — 4 times smaller with the exact kinematics 
than the approximate calculation. We have found that 
the approximate kinematics yields similar results when 
x g ~ 5.7 a; with the proportionality coefficient being the 
slowly varying function of Q 2 . 

In Fig. 9 we present the comparison of the kx and 
collincar factorizations for low values of Q 2 « 2 — 8 GcV 2 . 
We see that the two computations differ more in this re- 
gion. The lowest values are given by the kx factorization 
approach with the exact kinematics while the highest val- 
ues are given by the calculations with an approximate 
kinematics. At the lowest bin, Q 2 = 2 GeV 2 , the differ- 
ences seem to be smaller. This is due to the fact that 
in this region the calculation is dominated by the con- 
tribution of the quarks and the non-perturbative input 
which is the same in the kx and the collinear factoriza- 
tion formulae in our approach. The range in x has been 
extended down to x = 10~ 6 to cover the LHeC kinematic 
region. 

In the lowest Q 2 region the quark contribution and the 
non-perturbative gluonic component becomes dominant. 
This is why the differences between the exact and the 
approximate kinematics are starting to become smaller 
in this region as the whole Fx gets larger contribution 



from the quarks which are treated in the same way in 
both of these calculations. 

In the larger x region, below 10~ 3 , we observe that the 
kinematical effects are more significant. Interestingly, the 
gluon contribution originating from kx > k<j at x ~ 0.01 
is very small for the scales Q 2 up to about 10 GeV 2 and 
this region is completely dominated by the quark and 
the non-perturbative contribution. Therefore, using the 
approximate kinematics leads to large overestimation of 
the perturbative gluon contribution in this region. 



VI. CONCLUSIONS 

In this paper we have computed the longitudinal struc- 
ture function Fx of the proton within the kx factorization 
framework, using the unified BFKL/DGLAP resumma- 
tion scheme for the unintegrated gluon density. Since we 
have only fitted our parameters to the F 2 data, the calcu- 
lations for Fx are absolute predictions. The calculations 
are consistent with the experimental data from HERA 
collider. 

We have analyzed the impact of the exact kinematics in 
the kx factorization scheme. The exact gluon kinematics 
is very important for the phcnomenological description of 
the data on Fx- In particular, it leads to larger differences 
than changing from the kx factorization to the collinear 
factorization scheme, at least in the available kinematical 
region of the HERA data. 

We have also shown that the kx factorization scheme 
with exact kinematics includes both the collinear and 
the dipole limits. The first one is recovered by assuming 
the strong ordering in the gluon and quark momenta, 
together with the exact gluon kinematics. The latter one 
is recovered in the limit when the fraction of the gluon's 
longitudinal momentum is set to be equal to Bjorken x. 

The precision of the available HERA data for Fx does 
not allow to uniquely discriminate between different ap- 
proaches. It is possible however that by lowering the 
value of Q 2 or extending the kinematic regime to lower 
values of x (like in the proposed LHeC collider) one 
can explore the differences between the presented frame- 
works. 
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FIG. 5: Comparison between the collinear and the kr factorization calculations with exact gluon kinematics. The data are from 
the ZEUS experiment [2]. The light quarks u, d, s are treated as massless, the charm quark mass is set to be m c = 1.5 GeV. The 
solid (red) line denotes the calculation using the k factorization with exact kinematics, the black (dashed) line is calculation 
using the collinear factorization with the massive charm quark. 
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FIG. 6: Comparison between the collinear and the kr factorization calculations with exact gluon kinematics. The data are 
from HI experiment [1]. The light quarks u,d,s are treated as massless, the charm quark mass is set to be m c = 1.5 GeV. 
The solid-red line denotes the calculation using the kr factorization with exact kinematics, the black-dashed line is calculation 
using the collinear factorization with the massive charm quark. 
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FIG. 7: Comparison between the exact and the approximate (dipole like) kinematics in the kr factorization formula. The 
data are from the ZEUS experiment [2]. The light quarks u,d,s are treated as massless, the charm quark mass is set to be 
m c — 1.5 GeV. The solid (red) line denotes the calculation using the kr factorization with exact kinematics; the dotted ( blue) 
line is calculation using the kr factorization with the approximate kinematics. 
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FIG. 8: Comparison between the exact and the approximate (dipole like) kinematics in the fcr factorization formula. The data 
are from HI experiment [1]. The light quarks u, d, s are treated as massless, the charm quark mass is set to be m c — 1.5 GeV. 
The solid, red line denotes the calculation using the kr factorization with exact kinematics; the blue, dotted line is calculation 
using the fcr factorization with the approximate (dipole-like) kinematics. 
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FIG. 9: Comparison between the collinear and the kr factorization calculations for low values of Q 2 . The range in x covers the 
LHeC kinematic range. The light quarks u, d, s are treated as massless, the charm quark mass is set to be m c = 1.5GeV. The 
solid (red) line denotes the calculation using the kr factorization with exact kinematics; the dotted (blue) line is calculation 
using the kr factorization with approximate kinematics and the dashed (black) line is the calculation using the collinear 
factorization. 



